WIMP diffusion in the Solar System including solar WIMP-nucleon scattering 
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Dark matter in the form of Weakly Interacting Massive Particles (WIMPs) can be captured by 
the Sun and the Earth, sink to their cores, annihilate and produce neutrinos that can be searched for 
with neutrino telescopes. The calculation of the capture rates of WIMPs in the Sun and especially 
the Earth are affected by large uncertainties coming mainly from effects of the planets in the Solar 
System, reducing the capture rates by up to an order of magnitude (or even more in some cases). 
We show that the WIMPs captured by weak scatterings in the Sun also constitute an important 
bound WIMP population in the Solar System. Taking this population and its interplay with the 
population bound through gravitational diffusion into account cancel the planetary effects on the 
capture rates, and the capture essentially proceeds as if the Sun and the Earth were free in the 
galactic halo. The neutrino signals from the Sun and the Earth are thus significantly higher than 
claimed in the scenarios with reduced capture rates. 



INTRODUCTION 

The identity of dark matter is one of the major prob- 
lems in physics. If WIMPs (weakly interacting massive 
particles) exist they naturally produce the correct relic 
abundance, making them one of the best dark matter 
candidates. There are many ongoing searches hoping to 
observe the particle nature of dark matter (see [T]) for a 
review) . Searches for neutrinos from WIMP annihilations 
in the Sun and the Earth are receiving a lot of interest 
with the recent completion of the IceCube neutrino tele- 
scope and the recently proposed dark disc component of 
the Milky Way dark matter halo [SJ|3]. 

The idea to search for signals from dark matter accu- 
mulating in Solar System objects by WIMP-nucleon scat- 
terings, so called weak capture, goes back to [?]. Even if 
the main ideas were laid out already in the mid 1980s, 
consensus has still not been reached for the capture rate 
of WIMPs in the Sun and the Earth. The main compli- 
cation is that the Sun and the Earth are both part of the 
Solar System. 

For the Earth, things are complicated by the Earth 
being deep within the potential well of the Sun, accel- 
erating the galactic WIMPs reaching the Earth. The 
history of determining the capture rate has gone through 
a series of rather abrupt changes as our understanding of 
the problem has evolved (for a review, see [5])- For res- 
onances, i.e. when the WIMP mass matches the mass of 
an element in the Earth, the Earth can efficiently capture 
WIMPs directly from the galactic halo as in [B]. It was 
later realized that far from the resonances, i.e. for heavy 
WIMPs, the Earth's low escape velocity only allows it to 
capture WIMPs of velocities lower than the solar escape 



velocity at the Earth's distance r®, i.e. only WIMPs al- 
ready bound to the Solar System. It was however found 
in [7] that gravitational interactions with the planets will 
cause WIMPs to diffuse into the Solar System, thus form- 
ing a bound population of WIMPs. The gravitational 
interaction by the planets diffuse WIMPs between the 
population bound to the Solar System and the galactic 
WIMP population; given enough time this will result in 
an equilibrium configuration of WIMPs in the Solar Sys- 
tem, referred to as gravitational equilibrium. The grav- 
itational equilibrium configuration is easily found since 
the gravitational interactions causing the diffusion obey 
Liouville's theorem, assuring that the diffusion preserves 
the phase space density. Hence, any WIMP phase space 
over- or under-density in the Solar System will, given 
enough time, be mixed with the large galactic WIMP 
reservoir and recover the equilibrium configuration, i.e. 
the same phase space density of Solar System bound 
WIMPs as for the galactic WIMP halo [7 . One could 
also argue for this from detailed balance. Since the dif- 
fusion in and out of the Solar System is symmetric, we 
will at some point reach an equilibrium between parti- 
cles diffusing into and out of the Solar System, and this 
equilibrium will happen when the phase space densities 
are the same. 

It was later realized [8] that the bound WIMP popu- 
lation could be reduced via solar depletion, which occurs 
if WIMPs in the bound population are perturbed to end 
up in the solar loss cone, i.e. on orbits crossing the Sun, 
allowing them to be captured by the Sun. If this process 
is more efficient than the gravitational diffusion refilling 
the bound population, a reduction could occur. This was 
studied with numerical simulations in [5j and a reduction 
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of the capture rate for heavy WIMPs of up to an order 
of magnitude was found (shown as the magenta, dotted 
curve in Fig. [2j. 

We also mention a bound WIMP population proposed 
in [5] emerging from weakly captured WIMPs whose or- 
bits only pass the outskirts of the Sun. These WIMPs 
can be gravitationally disturbed by the planets onto or- 
bits no longer crossing the Sun. It was argued in [S] that 
this WIMP population will over the lifetime of the Sun 
build up a phase space density much higher than that of 
the galactic halo, making it important for WIMP cap- 
ture in the Earth. Such gravitational interactions must 
obey Liouville's theorem, making it impossible for the 
proposed gravitational diffusion to create a phase space 
density larger than that of the solar loss cone. The latter 
will be discussed in detail later, for now we just conclude 
that the high phase space density of |9J cannot exist. This 
population was also numerically investigated in [10], find- 
ing no support for a large enhancement. 

For the Sun, the capture rate calculations are simpler 
as the Sun gravitationally dominates the Solar System. 
However, the WIMPs in the solar loss cone can be dis- 
turbed gravitationally by the planets, and it was shown 
in [11 that the WIMPs reaching out to Jupiter will be 
thrown out by Jupiter before having time to scatter again 
in the Sun, introducing a Jupiter depletion. For WIMPs 
being heavy enough compared to the target nucleus, the 
energy loss in the scatter is small enough for the first 
bound orbit of all weakly captured WIMPs to reach ro_ 
and hence be perturbed before having time to scatter 
again in the Sun. The scattering in the Sun occurs 
either via spin-dependent scattering (mainly on hydro- 
gen) or spin-independent scattering (also on heavier ele- 
ments). In [TT] it was concluded that Jupiter depletion 
substantially reduces the solar capture rate for WIMPs 
heavier than a TeV, which scatter predominantly spin- 
dependently in the Sun. 

In this paper, we will show that Liouville's theorem is 
not only valid for gravitational diffusion but also holds 
approximately for weak capture. To a first approxima- 
tion the phase space density in the solar loss cone is then 
the same as that of the gravitationally captured WIMPs. 
Hence, weak capture should rather be viewed as another 
source of WIMP diffusion in the Solar System. For scat- 
tering cross sections that are not too low, the equilibrium 
time scale for the solar loss cone is shorter than the age of 
the Solar System. For a 100 GeV WIMP, equilibrium has 
been reached for spin-dependent scattering cross sections 
(on protons) ctsd ^ 10~ 44 cm 2 or spin-independent scat- 
tering cross sections (on protons) asi ^ 10 -46 cm 2 . For 
heavier WIMPs, slightly larger cross sections are needed 
[111 fT2"]. However, not all regions of bound phase space 
have reached gravitational equilibrium, giving a region of 
lower phase space density, a hole, as shown in Fig. 3 in 
[TJ. This so called Gould hole affects velocities (with re- 
spect to the Earth) between ~ 27 km/s and ~ 69 km/s. 



Weak capture will populate a vertical strip at ~ 30 km/s 
in this hole. 

The discussed solar and Jupiter depletions are both, in 
principle, gravitational diffusion between the two bound 
populations; as their phase space densities are approx- 
imately the same the effects, to a first approximation, 
cancel completely as the net flow between them vanishes. 
The studies of [TT] still found these effects important 
since they only looked at diffusion in one direction, ne- 
glecting the reverse process. In [5] they assumed the Sun 
to be optically thick to WIMPs, i.e. removed all WIMPs 
entering the solar loss cone, neglecting diffusion out of 
the solar loss cone. In the Jupiter depletion calculation 
of [TT] the possibility of an inverse process, i.e. diffusion 
into the solar loss cone, was not considered. In the subse- 
quent paper (T3] , diffusion of WIMPs by Jupiter into the 
solar loss cone was actually observed but not investigated 
further. 



LIOUVILLE'S THEOREM AND WIMPS IN THE 
SOLAR LOSS CONE 



For Liouville's theorem to be applicable, the force, F, 
encountered by the particles must be differentiable (i.e. 
"smooth") and must fulfill V p ■ F = where V p is the 
gradient operator in momentum space. Hence, Liouville's 
theorem is applicable to forces that do not depend on 
momentum. Obviously, weak scattering satisfies neither 
of these requirements. However, we will see that weak 
capture can be well approximated by a force fulfilling 
Liouville's theorem. 

WIMPs passing through the Sun typically scatter or 
1 times per passage. In such a scatter the WIMP loses 
the energy (over WIMP mass) AE uniformly distributed 
in the range 



< AE < 



AMm 
(M + m)' 



;Ek 



(1) 



where M and m are the masses of the WIMP and the 
target nucleus, respectively, and E^ is the kinetic energy 
over WIMP mass of the WIMP just before the scatter. 
The WIMP scatter probability in the solar passage de- 
pends very weakly on the WIMP momentum, but the 
energy loss in the scatter goes as AE oc Ek, as seen in 
Eq. 0. To fulfill Liouville's theorem the WIMP energy 
loss should rather be proportional to the time the WIMP 

— 1/2 

spends inside the Sun, i.e. AE oc E k . 

Since we ultimately want to investigate gravitational 
effects of the planets, we are only interested in WIMPs 
in the galactic halo or in the solar loss cone on orbits 
reaching out to the planets. The innermost planet in the 
Solar System is Mercury, at radius rj-j. A bound solar 
crossing WIMP on an orbit stretching out to at least 
will at radius r inside the Sun have kinetic energy in the 
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range 



tf (r) 1 




< E k < * (r), fc 



Mr) 
* (i? ) 



> 1, 



(2) 

where i? is the radius of the Sun and \P (r) is the ab- 
solute value of the solar gravitational potential at the 
distance from the Sun r. As R Q ~ 0.012r^, the WIMP 
kinetic energy inside the Sun is well approximated by 
Ek — especially for orbits going further out to 

the more massive and more interesting planets. Further- 
more, the typical WIMP velocities in the galactic halo 
are significantly lower than the escape velocity inside the 
Sun and unless we are on a resonance, only the low ve- 
locity WIMPs of the halo are captureable by the Sun, 
making Ek ~ ^©(V) a good approximation also here. 
Since all WIMPs of interest to us have, essentially, the 
same velocity inside the Sun, the momentum dependence 
in the scatter is not important for the WIMPs of interest 
except, possibly, if we are close to a resonance. 

For an outside observer, who does not notice the 
small momentum dependence, the random energy loss of 
WIMPs passing the Sun could be equally well-attributed 
to a time dependent deepening potential well inside the 
Sun, i.e. a differentiable force fulfilling Liouville's theo- 
rem. Hence we conclude that Liouville's theorem can be 
applied to high precision for the weakly captured WIMPs 
that are subject to gravitational interaction with the 
planets. 

In the language of detailed balance, we can see this 
approximate relationship from the fact that the rate at 
which WIMPs scatter into the solar loss cone is approxi- 
mately the same as the rate at which they leave via fur- 
ther scatterings in the Sun (as the velocity dependence 
is so weak). This means that the phase space density of 
WIMPs in the solar loss cone will be approximately the 
same as in the galactic halo. As we move close to the 
Sun (within the radius of Mercury), the approximation 
of equal rates worsens (Liouville's theorem is less applica- 
ble), and the density will be different. We are, however, 
not concerned with those orbits. 

We will now verify the above reasoning by calculating 
the density in the solar loss cone using Liouville's theo- 
rem and compare with the numerical results of |12) . We 
assume a Maxwell-Boltzmann distribution for the galac- 
tic halo as in |12| but this will not be important for our 
conclusions. The same arguments will hold also for other 
velocity distributions. The relevant galactic phase space 
density is then, viewed by an observer (i.e. the Sun) mov- 
ing with the rotational velocity u = y/2/3v, given by 



16. 



F(u) 



n w - 



T 3/2 



3/2 



exp 



3m 2 
ZD* 



sinh(\/6u/w) 
(3) 



Here u = |u| and n\y are the local, galactic WIMP ve- 
locity and number density far from the Sun, respectively; 
v = 270 km/s is the three-dimensional velocity disper- 
sion. F(u) maximizes as u — > and is approximately 
constant for low velocities, i.e. for u -C v. Since only 
the very low velocity WIMPs in the galactic halo are 
prone to planetary gravitational diffusion, the gravita- 
tional equilibrium phase space density for WIMPs in the 
Solar System is F(0) [TJ, which is also the equilibrium 
configuration with maximal density. The maximal phase 
space density in the solar loss cone is then also F(0) and 
is achieved if weak capture only probes very low it, i.e. for 
M >mas seen in Eq. (JT|). 

Knowing the phase space density, the spatial density in 
the solar loss cone is simply the integral over the allowed 
angular momentum for velocities below the local escape 
velocity v e (r) with r being the distance from the solar 
centre. For a bound WIMP's orbit to cross the Sun, 
the angular momentum over WIMP mass J must fulfill 
J < Jmax = wRq, where w is the local velocity at i? . 
As for Ek in Eq. J max is essentially the same for 
all WIMPs of interest, giving J max ~ v /2GA/ i? Q . One 
can show that for an isotropic velocity distribution the 
probability for a bound WIMP at radius r with velocity 
u to be on a solar-crossing orbit, i.e. to fulfill J < J max , is 
P Q {u, r) = 1 - y/l - (J max /ru) 2 for ru > J max and P = 
1 otherwise. The maximal number density of WIMPs in 
the solar loss cone from Liouviullc's theorem is then 



4iru 2 F(0)P & {u,r) du 



6\/3 1 / GM(r 



ZypK V 3 \ R(. 



3/2 



Rr. 



-5/2 



(4) 



where we have made a Taylor expansion in i? jr which 
gives an error < 0.3% for r > ry. This result is shown 

in Fig. [I] (green, dashed curve) and compared with the 
detailed numerical calculation in |12| . They fit very well 
when the assumption M ^> m is fulfilled (which is al- 
ways the case for spin-dependent scattering). The nu- 
merical simulations in |12j makes no assumptions on Li- 
ouville's theorem and consists of a Monte Carlo study 
following WIMP orbits to build up the bound popula- 
tion. That the agreement with our theoretical calcula- 
tions is so good verifies that Liouville's theorem indeed 
applies here. Note that this figure is independent of the 
magnitude of the scattering cross sections. (It does de- 
pend on if capture is dominated by spin-independent or 
spin-dependent scatterings though.) See |12j for more 
details. 

If M ~ m, the Sun also captures galactic halo WIMPs 
of higher velocities, which renders F(u) ~ F(0) (i.e. u <C 
v) a bad approximation. Taking this into account reduces 
the normalization of the solar loss cone density with the 
factor F(u)/F(0) averaged over the galactic velocities u 
of the captured WIMPs (which reach the planets). We 
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FIG. 1. The WIMP line number density n(r) = <t-Kr 2 n%{r) 
in the solar loss cone for different WIMP masses with spin- 
dependent scatter dominating in the Sun (left figure) and 
spin- independent scatter dominating (right figure). The 
green, long-dashed lines are the maximal densities by Liou- 
ville's theorem, i.e. for M 3> m, as in Eq. 0. The blue, 
dotted lines show the density by Liouville's theorem but tak- 
ing into account that F(u) is not flat and is hence valid also 
for M ~ m. These are compared with the numerical results 
from SE 2009 |12j shown with red, solid lines. 



have performed this calculation numerically and the re- 
sult is shown in Fig. [I] (blue, dotted line), giving a better 
fit to the numerical results for M ~ m. We see that 
our prediction fits very well with the numerical results 
with a slightly worse fit close to the iron resonance, i.e. 
for spin-independent scattering with M — 100 GeV. We 
believe that this is due to Liouville's theorem not being 
a good approximation for the first scatter so close to the 
resonance since then all halo velocities of WIMPs bound 
to the Milky Way are prone to weak capture, making 
Ek ~ ^0 (r) of Eq. ^ a less good approximation. Spin- 
independent capture of light WIMPs in the Sun occur 
mainly on helium, carbon, nitrogen, oxygen and iron. As 
the masses vary, there will always be some significant 
capture not close to a resonance, the error in the blue 
curve in Fig. [T] is thus not expected to be much larger for 
lighter WIMPs. Also, the bound population is not very 
important close to the resonances (see Fig. [2) , hence the 
discrepancy is not very important. For smaller radii than 
the planets, the agreement in Fig. [I] is worse, but this is 
expected as several of our approximations break down in 
the limit r ~ R@; this is of no concern to us as we are 
only interested in the WIMP density at r > rw. 



EARTH'S WIMP CAPTURE RATE AND SOLAR 
DEPLETION 

To make a conservative estimate on how much solar 
depletion can affect the Earth's capture rate, we assume 
the worst case scenario, i.e. that the refilling of the solar 
loss cone and the diffusion between the two populations 
are so efficient that the phase space density of all bound 
WIMPs is reduced to that of the solar loss cone. The 
capture rate for this conservative case is shown in Fig. [2] 
concluding that even this maximal effect is indeed very 
small, at most a reduction by ~ 35% for specific WIMP 
masses. The reduction is even smaller if spin-dependent 
scatter dominates in the Sun. The reason for the small 
reduction is that, for low enough WIMP masses to signif- 
icantly reduce the solar loss cone density, the Earth effi- 
ciently captures WIMPs directly from the galactic halo. 
As equilibrium between capture and annihilation has typ- 
ically not been reached in the Earth, the annihilation rate 
goes as the square of the capture rate (see e.g. [S]), mak- 
ing the boost of the annihilation rate even larger com- 
pared to [5j. 

In Fig. [2] we also show the capture rate for the 'stan- 
dard' approximation |7J of taking the velocity distribu- 
tion as if the Earth were in free space (i.e. at the location 
of the Solar System, but not affected by the gravitational 
attraction by the Sun) . Also shown are the capture rates 
directly from the galactic halo, i.e. the unbound popula- 
tion, as well as from the population bound to the Solar 
System with the Gould hole in Fig. 3 in (7] . From Fig. [2] 
we see that the capture rate is not very sensitive to so- 
lar depletion or the Gould hole, concluding that the free 
space approximation is actually quite good. 

For clarity, we also show the velocity distributions for 
different WIMP populations that go into the capture rate 
calculations. We can define the one-dimensional velocity 
distribution f(u) as f(u) = ^-F(u) with F(u) from 
Eq. (|3|. In the capture rate calculations, f(u) enters 
as f(u)/u (see e.g. [FJ [7] for details), which is shown in 
Fig. [3] We can clearly see that in the absence of the 
Gould hole, the bound plus unbound population almost 
perfectly matches the free space velocity distribution, as 
argued above. The Gould hole reduces the velocity distri- 
butions between ~ 27 km/s and ~ 69 km/s. The bound 
population is also affected by weak capture with the pos- 
sible reduction discussed above. 

We have assumed a Maxwell-Boltzmann velocity dis- 
tribution, but our method applies to any velocity dis- 
tribution. If the velocity distribution is very peaked at 
low velocities, the reduction in the solar loss cone could 
be larger, i.e. the difference between the green, dashed 
curves and blue, dotted curves in Fig. [T] could be larger. 
Inversely, if the velocity distribution is peaked at high ve- 
locities (for a stream e.g.), we could get an enhancement 
instead of a reduction from weak capture. 
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FIG. 2. The capture rate in the Earth for different phase space densities. We show the Gaussian (Maxwell-Boltzmann) at the 
Sun's location but for the Earth as in free space, the "best" curve from [Sj, the unbound population, the bound population with 
Gould's hole empty, the total with the hole and the total with the maximal reduction of the bound densities for spin-dependent 
(SD) and spin-independent (SI) scattering respectively dominating in the Sun. In the left figure, we show the capture rates and 
in the right the ratios with respect to the free space Gaussian. The cross section chosen is just for reference and the capture 
rate scales linearly with it. 



PLANETARY EFFECTS ON THE SUN'S 
CAPTURE RATE 

Jupiter is with its mass of 318 Earth masses the planet 
in the Solar System with by-far the largest gravitational 
impact, as observed in [TT] in removing WIMPs from 
the solar loss cone. The other consequence of this is that 
Jupiter gravitationally diffuses WIMPs efficiently enough 
to have achieved gravitational equilibrium for all parts of 
phase space representing bound orbits crossing 0[7]. 
Jupiter depletion is only important for M 3> m but the 
solar loss cone phase space density is then the same as 
that of gravitational equilibrium. Hence, Jupiter by Li- 
ouville's theorem does not alter the solar loss cone density 
at since WIMPs are equally diffused into as out of the 
solar loss cone, leaving no net effect on the Sun's WIMP 
capture rate. 

Another way to see this is to calculate the solar capture 
rate from instead of from r — > oo. As gravitational 
equilibrium has been reached at r\, the total phase space 
density (bound plus unbound WIMPs) is the same at 
as in the galactic halo, i.e. F(rq,,u) = F(oo, u) = ^(u). 
The solar WIMP capture rate from is then the rate of 
WIMPs scattering in the Sun to orbits with r max < ro_, 
analogous to the standard scenario where the WIMPs 
from the galactic halo scatter to an orbit with r max < oo. 
The deviation of the capture rate from compared to 
the standard scenario can be incorporated by assigning 



an "effective" escape velocity v e defined as the minimal 
velocity required at the location of the scatter to reach 
. Comparing v e with the usual escape velocity v e gives 

- = \l^¥> ® 

where k > 1 depends on the location of the scatter and is 
defined in Eq. Q. As ~ 1119i?0 we get that v e ~ v e , 
concluding that the Sun's capture rate is essentially un- 
affected by the gravitational presence of Jupiter. This 
presents an alternative way of reaching the same conclu- 
sion as we did before, i.e. that the influence of Jupiter 
alone is not important for the solar WIMP capture rate. 

In this paper we have assumed the phase space den- 
sity at rq^_ to really be F(0), as in [7]. This is sup- 
ported by the numerical work in [£]. In jTHI, a slightly 
lower bound phase space density is found for some ve- 
locities. We argue that this is probably due to having 
neglected regions of longer equilibrium time scales (ret- 
rograde WIMPs with high angular momentum) as they 
are below the resolution limit of the simulations. It is 
argued in [13J that these WIMPs are not important by 
referring to Fig. 3 in |T3] , but also there the resolution of 
the relevant WIMPs is too low to draw any conclusions 
on this population. 
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a capture rate somewhere in between the conservative 
and the unreduced (total) capture rate (i.e. between the 
black curves and the dark green, dot-dot-dashed curve in 
Fig. [2|, depending on the relative efficiency of weak and 
gravitational diffusion. 

The only way to achieve a higher phase space density 
of bound WIMPs in the Solar System is to either have 
a mechanism that violates Liouville's theorem, to have a 
higher WIMP phase space density in the galactic halo, 
or if the Solar System in some way was born with a dark 
matter overdensity (of which a small portion, e.g. in the 
Gould hole, could survive until today). The second case 
is achieved by the galactic dark matter having a disc 
structure which co-rotates with the stars, increasing the 
phase space density of the important low velocity WIMPs 

[HE]. 
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FIG. 3. The low velocity part of the velocity distributions (di- 
vided by velocity) of the various unbound and bound popula- 
tions assuming a standard Maxwell-Boltzmann velocity distri- 
bution in the halo. Without inclusion of the Gould hole, the 
bound + unbound population almost perfectly matches the 
free space approximation. The Gould hole (i.e. taking into 
account that part of the phase space for the bound popula- 
tion is not filled), affect intermediate velocities and therefore 
has no dramatic effect on the capture rates shown in Fig. [2] 



SUMMARY AND DISCUSSION 

We have shown that one can to high precision also 
use Liouville's theorem for weakly captured WIMPs, not 
just for the gravitationally captured WIMPs as previ- 
ously believed. The solar weak capture process is hence 
to be viewed as an extra source of WIMP diffusion in the 
Solar System. 

Both the effects of solar and Jupiter depletion have 
in the literature been found to be the most pronounced 
for heavy WIMPs, i.e. when the phase space density of 
the solar loss cone is the same as that of gravitational 
equilibrium. As both these effects are really gravita- 
tional diffusion between the solar loss cone and the pop- 
ulation bound through gravitational diffusion, the effects 
are by Liouville's theorem cancelled by diffusion in the 
reverse direction. This results in the overall conclusion 
that WIMP capture in the Sun and the Earth can gen- 
erally be treated as if they are both free in the galac- 
tic halo, i.e. returning to the conclusion of [7]. Only 
the Earth's capture rate can deviate from this with, at 
most, 30-35% where the most important factor is the na- 
ture of weak capture in the Sun. Note that this is under 
the conservative assumption of weak capture being much 
faster than gravitational diffusion. In reality we expect 
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